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The superconducting upper critical field Hc 2 {T) of a two dimensional BCS superconductor is calcu- 
lated in the vicinity of a van-Hove singularity. The zero temperature upper critical field is strongly 
enhanced at weak coupling when the Fermi contour coincides with van-Hove points, scaling as 
Hc2{0) oc compared to the usual result Hc2{0) oc . The result can be interpreted in terms of 
the non- Fermi liquid decay of normal state pair correlations in the vicinity of a van-Hove point. 
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The possibility of a strong enhancement of the super- 
conducting transition temperature in the vicinity of van- 
Hove singularities (vHs) oL electronic systems has been 
raised on several occasionsE’cl At a vHs, the Fermi veloc- 
ity vanishes at points on the Fermi contour and in 2D, 
the density of states diverges logarithmically. While it is 
tempting to apply this simple model of U enhancement 
to realizations of «tiasi-2D systems such as Copper-Oxide 
superconductors,® it is often pointed out that three- 
dimensionality and scattering effects have the effect of 
smearing out the vHs. Also, the electron kinetic energy 
favors low density of states at the Fermi level and struc- 
tural distortions may be expected to occur in order to 
push the vHs below the Fermi level.B 

Despite these limitations, the 2D vHs provides a sim- 
ple and well-understood example of a normal Fermi sys- 
tem which possesses anomalous normal state correla- 
tions and therefore deserves to be understood in de- 
tail. The origin of non-Fermi liquid effects in this case 
is the absence of a well-defined Fermi velocity scale, 
which is in turn due to the presence of an underly- 
ing lattice and the absence of Galilean invariance. It 
has recently been shown that non-Fermi liquid behavior 
in the normal state leads to deviations from the usual 
quasi-parabolic mean- field Hc 2 (T) curve of Fermi Iku 
uid BGS (FL-BGS) superconductors found by Gorkov.® 
Interestingly, deviations from the parabolic jshape have 
been observed in layered superconductors,u includiBg 
both Gopper-Oxide^ andi-(BEDT-TTF) 2 Gu(NGS) 2 .E^ 
As emphasized previously^ the mean field upper critical 
field is a normal state property which probes both the 
spatial (via the magnetic length) and energy (tempera- 
ture) dependence of pair correlations. It may therefore 
prove an effective probe of non-Fermi liquid normal ef- 
fects. 

In the present paper we show that anomalous corre- 
lations in the vicinity of a 2D vHs lead to spectacu- 
lar deviations from the FL-BC^ result. Following the 
semi-classical Gorkov approach,Ej’E3 the superconducting 
transition corresponds to the appearance of non-trivial 



solutions of the equation 

A{r)=g (I) 

where j3 is the inverse temperature and g the pair- 
ing interaction. The magnetic field effect is ip-i 
eluded in a semi-classical form as a phase termlJ 
exp(z[A(r) -b A(r')](r' — r)). This approximation ne- 
glects Landau level quantitization which is expected to 
become important in the high field limit. The latter be- 
havior has not yet been observed experimentally within 
the range of fields presently accessible and the semi- 
classical approximation will be used in this paper. K^ir) 
is the normal state pair propagator in the absence of an 
external magnetic field: 

K/3{r) = f3~^'^g{r,ujn)g{r,-ujn)- ( 2 ) 

CJ-n 

Here g(r,Wn) is the real space single particle green func- 
tion at Matsubara frequency The highest eigenvalue 
of Eq. © determines the upper critical field. We work 
in the Landau gauge A — (0, Kxy 0). Making use of the 
degeneracy of the gap function,Eil one can write 

A{x) = g J dx'Kis[x' — x,—H{x + x')]A{x'), (3) 

where A(a;) is the gap function integrated over y and 
Ki^{x, ky) is the Fourier transform of Kj^^x, y). The latter 
can be expressed as 

Kij{x,ky) = 2 J dco J dpy tanh(/3w/2) x 

A{x,Py + ky,ui)B{x,-py,-uj), (4) 

where A{x,ky,oj) = ^ / dfca, exp(ifca,a:)Img(fc, w) and 
B{x,ky,uj) = -^ J dkx exp{ikxX)Keg{k,uj). 

When the Fermi energy Ep is close enough to the van- 
Hove singularity one can express E{k) — E'"^ = qxQy 
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FIG. 1. Constant energy contours of the quasiparticle dis- 
persion near a 2D van Hove point E{k) — = <lx<ly 

The values of q are restricted to a region around the saddle 
point by the cutoffs Po and No- 



where q — k — and the values of the momentum are 
restricted to a small region around the saddle point 
by a cutoff. For simplicity, we assume that 2k""^ is a 
reciprocal lattice vector, but the results are independent 
of this assumption. We now have 



A(x,qy,u) = - 



1 x(u: + 5) 

e e 
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B{x,qy,uj) = isgn{ — )A{x,qy,u;), 

qy 



( 5 ) 



with S = Ep — The pair propagator is 
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( 6 ) 



We have introduced a cutoff Po for the y component of 
the momentum. Note that doping away from the vHs 
introduces a oscillating term in the integral and Kpix^ ky) 
decreases. The effects of the temperature and the doping 
in the integral are similar as they provide a cutoff for 
large x. One expects Hc 2 (T) to be rather insensitive to 
T for T < S, that is, i?c 2 remains approximately constant 
to a temperature of the order of <5. For T > S, we have the 
opposite, that is, the i?c 2 curves are fairly independent 
of the filling. At zero temperature the asymptotic decay 
of Eq. (^) in real space is AT ~ This is a slow 

decay relative to the 2D Fermi Liquid (circular Fermi 
surface) result K ~ {x“^ +y‘^)~^. 

Firstly, we evaluate the critical temperature at zero 
field when the vHs lies at the Fermi level. We-.need the 
q = 0 pair propagator (or pair susceptibility) :llJ 
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FIG. 2. Variation of the Hc2 curve as a function of S, for 
Pojxo = 20 and moderate coupling g — 0.2, obtained by nu- 
merical lanczos solution of Eq. (3). Notice that Hco has a 
much stronger dependence on 5 that Tco- Upward curvature 
is observed in all curves down to a temperature of the order of 
5. For higher 5, there is a small increase of JTc 2 with temper- 
ature as the temperature averaged density of states increases. 



Here we have, introduced a short distance cutoff Xg in the 
x-direction.E3 At temperatures low compared with the 
bandwidth XoT^/Ppo 1 , the main contribution to the 
integral results from the small x region and /C/ 3 ( 0 , 0 ) ~ 
1/2 log^ [xoTt/ I dpo]- The zero field critical temperature is 
determined by the relation 1 = gK.p^ ( 0 , 0 ) and one re- 
covers the standard result Tc = n~^x~^Po exp(— -^ 2 / 5 ) 
compared to the weak coupling FL-BCS result T/ oc 
exp(— 1 /^, where g fx g. 

The pair propagator [Eq. ®] can be evaluated in 
closed form at zero temperature and zero S. In this case, 
the gap equation [Eq. (||)[ simplifies to 



A{x) =g dx' 



2\x' - 



■log 



E[\x' + x\ 



A(x'), (7) 



with \x' — x\ > Xo- A lower bound for Elco can be ob- 
tained by the variational method. A physically reason- 
able form for the gap function is A(x) = l/v^B(a— |x|), 
where a is a length to be determined. The maxi- 
mum eigenvalue is obtained for a = e^XoPo/‘^H which 
leads to 1/g > 1/4 log^ {e^Po/2.XoHco) ■ Thus iLco 
[e^Po/2xo)exp[—2g~^). This implies 



^o^co f ^o-^co 

Po \ Po 

This result should be contrasted with the usual FL-BCS 
scaling Elco ~ T[?o- By expanding the kernel [Eq. ®] 
around the zero temperature critical point, one can ob- 
tain the low temperature behavior of i 7 c 2 and behaves as 
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Hc-Hco ~ -T\og{Hco/T)/log{xoHco/Po)- Hc{T) does 
not saturate as T goes to zero. 

For temperatures close to the zero field critical temper- 
ature Tco, one can follow the Gorkov microscopic deriva- 
tion of the Ginzburg-Landau equation to obtain 

aA{x) + b{d^ — H^x^)A{x) = 0, (9) 



where 






d^rKi3{r 






6 = i J d‘^rK/3{r)x‘^ ^ i j d'^rKf3{r)y'^. 



Here, we have used the fact that K/ 3 (r) is symmetric in 
the transformations ^ y. Setting po ^ then bx ^ by. 
For our van-Hove model with (5 = 0, we have 
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(10) 



Note that for the 2D isotropic BCS superconductor, we 
would have a ^ p{Ep){Tc — T)/Tc- The vHs at the Fermi 
level is averaged out at a finite temperature Tco and so, 
in the previous expressions, the density of states is re- 
placed by an effective density of states \og[T:XoTco/po\ 
giving a weak enhancement of the slope. Eqs. and 
show that there is an anomalous enhancement and 
consequently upward curvature of the upper critical field 
and that this upward curvature is strongly enhanced at 
weak coupling (low Tc). Thus in contrast to FL-BGS, 



normalized plots of (^) not fall on to 



a univer- 



sal curve for different couplings. 

To extend the above analysis to all temperatures and 
confirm the variational analysis requires a numerical ap- 
proach. Discretizing in the x direction, Xn = ne where n 
are integers, the integral gap equation becomes 



Ai — g'^^KijAj. ( 11 ) 

b 



While the system is now discrete it remains infinite di- 
mensional. However, the system can be truncated for fi- 
nite H because we need only consider solutions which 
are localized about the origin on a scale of the magnetic 
length. 

The highest eigenvalue and corresponding gap func- 
tion of Eq. were obtained by the Lanczos method. 
The results are displayed in Figs. || and In Fig. ||, 
we have Hc 2 curves for different values of 5 and fixed 
coupling g. Hc 2 saturates when T < 5. When i5 = 0, 
that is, when the Fermi energy coincides with the vHs, 




FIG. 3. Gap functions A(a;) for two points of the Hc 2 curve 
with poixo = 20, g = .2 and S — 0. One of the points is near 
the top of the Nc 2 curve and the other near the bottom. Al- 
though the range of the gap functions is roughly the same 
(of the order of the magnetic length in each case), the func- 
tional forms of the two curves are different. For low Hc/Hco 
we have the conventional behavior A(®) ~ e , while for 
Hc/Hco ~ 1, the gap function has an anomalously slow decay. 



there is no saturation and the Hc 2 curve shows an up- 
ward curvature through the complete temperature range 
in complete agreement with the variational calculation. 

A qualitative understanding of these results is ob- 
tained as follows. In FL-BGS, the upper critical field 
is given by 1 = gp{Ep)\og{H -\- T^). In the vHs case, 
p{E) is logarithmic divergent. Finite magnetic field or 
temperature smear out the density of states and p{Ep) 
should be replaced by log(H -I- T) at zero 6. Thus 
1 ~ g\og{H -\- T) \og{H -\- r^). This leads to Hco ~ Tf^ 
as found above. 

In conclusion, we have studied the suppression of the 
mean field superconducting instability of a clean weak 
coupling BCS superconductor in finite magnetic field at 
a 2D vHs. The upper critical field is strongly enhanced 
relative to Tc, this enhancement is described by the re- 
lation i?c2(0) ^ t/^. The upper critical field falls lin- 
early with temperature near T = 0. These effects disap- 
pear rapidly when the Fermi level is tuned away from the 
Fermi level; more precisely, they are absent when Tc < 5 
defined above. The 2D vHs provides a simple example 
of a system where non-Fermi liquid normal state corre- 
lations show up strongly in the temperature dependence 
of the superconducting upper critical field. 
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There is some freedom in the choice of the phase term. An 
equivalent expression would be exp{i[2A(r)](r' — r)} and it 
would lead to the same results. 

An equivalent expression is 

/C(0,0)= [ dpy— log 

Jo Py 

where the integrand is the 2D isotropic uniform pair prop- 
agator with Fermi velocity vf = Py- 

It can be shown that Hc 2 {T) depends only on the ratio 
PojXo- 
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